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Abstract 

We consider the 2D critical Ising model on a strip with fixed boundary conditions. It is 
shown that for a suitable reparametrization of the known Boltzmann weights the trans- 
fer matrix becomes a polynomial in the variable esc (4m), being u the spectral parameter. 
The coefficients of this polynomial are decomposed on the fixed boundaries Temperley-Lieb 
Algebra by introducing a lattice version of the Local Integrals of Motion. 

1 Introduction 

The objective of this paper is to extend the results of the previous work on Critical Dense 
Polymers [7J and Ising on a cylinder [5] to the case of Ising on a strip, in order to put the 
previous work on a more solid ground and to show that once more the same approach is 
possible. 

The approach goes all the way back to [I] , where the integrable structure of conformal field 
theory (CFT) was analyzed from the point of view of the existence of infinitely many higher 
spin conserved quantities in mutual involution. Inspired by this work we seek and succeed 
to obtain the tower of commuting conserved quantities directly on the lattice. 
The involutive charges of CFT are defined as the maximal Abelian subalgebra of the envelop- 
ing algebra UVir of the Virasoro algebra, the lattice analogue of this structure is obtained 
by replacing the Virasoro algebra with the Temperley-Lieb algebra, as it is well known [9]. 
There are of course different types of Temperley-Lieb (TL) algebras, in particular it makes 
a lot of difference to work with open boundaries or periodic boundaries [5] . In [7] we worked 
with an open boundary TL algebra, and obtained closed form expressions for many of the 
involutive charges, an ansatz closed form for the generic expressions was also provided in 
which some of the coefficients are still unfortunately not known. The heart of the difficulty, 
as pointed out long ago also in [H] lies in the presence of boundary terms which manifest 
themselves through the appearence of boundary tangles. This should be the case also for the 
Ising model on a strip, and infact we succeed in finding the same type of contribution from 
the boundary as we did for Polymers on a strip [2] [7J. We also succeed in finding a larger 
number of conserved charges on the lattice than we did for Polymers on the strip. 
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2 Generalities 



It is well known that the Ising model can be characterized in terms of Boltzmann weights, 
which can be derived from an elementary lattice hamiltonian [TJ [3] . Naturally such weights 
can be suitably rescaled by multiplying for trigonometric functions, in particular we shall 
employ here the following parametrization: 



W R I 

a 



u 



tan 



a 



) S a ,p + 8v-pj 8<t,-t + ( cot (t - uj 8a, P + 8 a - p j 8 a>T (2.1) 



cot(n)5 CTiP + 6 a -p^6 a> r + (tan(u)5 CTiP + 6 a - p j8 a - T (2.2) 



satisfying 



W R i 

a 



7T 

4 



u 



Wr, T 



a 



u 



(2.3) 



A transfer matrix can then be defined in different ways and for different boundary conditions, 
in particular if we choose to work on the strip we can build the transfer matrix in a boundary 
condition dependent way ,therefore we introduce the open boundaries transfer matrix: 
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Pi = l 

PL+2 = b 

and similar conditions for a. The transfer matrix has thus size 2 L x 2 L . 
The partition function is then defined as: 

Z = TrT A/ 



(2.4) 

(2.5) 
(2.6) 



(2.7) 



on general grounds the partition function is expected to behave asymptotically as a linear 
combination of CFT characters, whose modular parameter is related to the exponential of 
some trigonometric function of the spectral parameter. We will describe this behaviour of 
the partition function in a later section. 

It is also well known that the Ising model and in generally Spin-chains, Potts models, RSOS 
models, Vertex models [5] are directly related to representations of the Temperley-Lieb (TL) 
algebra. We will extensively investigate this connection for this special case, thus for our 
needs it is necessary to introduce the TL algebra, whose generators are defined by the 
following matrix elements: 



L+2 



{e 2 i)a,p = V28 aij(r . +1 Yl 8, 



(2.1 



fc=i 
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which satisfy the relations: 

,2 



= v 7 ^ (2.10) 



e^iie* = e, (2.11) 

e^e,- = ejei , |z - j\ > 2 (2.12) 

notice that there are 2L + 2 generators, due to the fact that the spin streaks are extended 
on the right and on the left boundary. 

3 Inversion Identity and Eigenvalues 

In this section we will solve the inversion identity for the transfer matrix thus obtaining the 
full spectrum of its eigenvalues. The approach here goes all the way back to [3], and more 
recently to pp. 

The inversion identity takes the following form: 

D(x)D(-x) = (-l) L+1 2 2L+4 t/ 2L+2 (x)l (3.1) 

where x = esc (4m) and U n (x) are the second kind Chebyshev polynomials. The Ui n can be 
factorized as: 

^ + ^-n>-sin*(<^)) (3.2, 



k=l 

L+1 



B(x)B(-x) = (_l)^+ 6 J] - sin 2 (^T^)) (3-3) 

k=i 

we then define the rescaled transfer matrix 

D(x) = 2 2L+3 x L+1 T(x) (3.4) 

hence the eigenvalues become 

k=l 

where \i = (p,i, ^■■■I^l, fJ>L+i) is such that /ij = ±1 one can alternatively introduce the 
partition V which is made up by the k such that ^ — — 1 and write 

x + sin(tfc) tt / sin(i 

kev 

where 



keV K ' k=l 



(2k-l)n . . 

The solutions of the inversion identity however, are too many, namely 2 L+1 vs the size of 
the transfer matrix which is 2 L . If we call p the length of the partition V we have that for 
b = +1 we must consider p as even, whereas for b = — 1 we must take p to be odd. This 
completes the discussion of the selection rules. 



3 



4 Temperley Lieb Decomposition 

In this section we want to decompose the transfer matrix in terms of the Temperley-Lieb 
algebra as: 

D(x) = 2 2L+3 ~ k ^x L+1 ~ k (4.1) 

rv • 

k=0 



and for the rescaled transfer matrix 

L+l 

k\(2x) k 



T (x) = Y -5* (4.2) 



k=0 



and taking the logarithm and expanding as x — > oo we define what we shall call the lattice 
integrals of motion A 2n _i 



k=l v ' 

being 

D fe = 5 fc (A 1 ,...,A fc ) (4.4) 
where the are the Bell polynomials, which are defined as 

EDO < 1° V — -V But^ , . 

fc=0 

the first A n that can be obtained from direct decomposition of the transfer matrix for low 
sizes read: 

A 2 = -(2L+1)1 (4.6) 

A 4 = -6(6L+1)1 (4.7) 

A 6 = -60(40L -4)1 (4.8) 

A 8 = -(352800L - 115920)1 (4.9) 

Aio = -(91445760L - 48625920)1 (4.10) 

A12 = -(36883123200L - 26424921600)1 (4.11) 

A M = -(21371135385600L - 18955051315200)1 (4.12) 

2L+2 

Ai = V2 ^ e, - (2L + 2)1 (4.13) 

i=l 

A 3 = H 3 + 6Ai - 2Bi + 41 (4.14) 

A 5 = 6H 5 + 60H 3 - 6[Ai, [Ai, Bi]] + 240Ai - 96Bi - 48B 2 + 2881 (4.15) 

A 7 = 90H 7 + 1260H 5 + 7560H 3 - 90[A 3 , [A l5 B^] - 1080[A 1; [A 1? B^] - 180[A 1; [A l5 B 2 ]]+ 
+ 25200Ax - 9360B! - 9360B 2 - 2160B 3 + 417601 

(4.16) 
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A 9 = 2520H 9 + 45360H 7 + 362880H 5 + 1693440H 3 - 420[A 5 , [A 1; B^] - 25200[A 3 , [A l7 Bi]] + 
- 5040[A 3 , [A 1; B 2 ]] - 342720[Ai, [Ai, Bi]] - 100800[A 1; [A 1; B 2 ]] - 20160[A 1; [A 1; B 3 ]] 
+ 5080320Ai - 1532160Bi - 2499840B 2 - 1048320B 3 - 161280B 4 + 104832001 

(4.17) 

An = 113400Hn + 2494800H 9 + 24948000H 7 + 149688000H 5 + 598752000H 3 - 1260[A 7 , [Ai, B x ]] 

- 189000[A 5 , [A l7 Bi]] - 37800[A 5 , [A 1; B 2 ]] - 9979200[A 3 , [A 1; BJ] - 3628800[A 3 , [A b B 2 ]] 

- 907200[A 3 , [Ai, B 3 ]] - 167832000[Ai, [Ai, Bi]] - 62596800[Ai, [Ai, B 2 ]]+ 

- 17236800[Ai, [Ai,B 3 ]] - 1814400[Ai, [Ai,B 4 ]] + 1676505600Ai 

- 355622400B! - 954374400B 2 - 555206400B 3 - 156038400B 4 - 18144000B 5 + 40787712001 

(4.18) 

A 13 = 7484400H 13 + 194594400H n + 2335132800H 9 + 17124307200H 7 + 85621536000H 5 + 
+ 308237529600H 3 - 2970[A 9 , [Ai, B x ]] - 831600[A 7 , [A u Bi]] - 166320[A 7 , [Ai, B 2 ]] + 

- 109771200[A 5 , [Ai, Bi]] - 39916800[A 5 , [Ai, B 2 ]] - 9979200[A 5 , [Ai, B 3 ]] + 

- 5688144000[A 3 , [A b Bi]] - 2634508800[A 3 , [A b B 2 ]] - 1017878400[A 3 , [A x , B 3 ]] + 

- 119750400[A 3 , [Ai,B 4 ]] - 114960384000[Ai, [Ai.Bj] - 49456915200[Ai, [Ai,B 2 ]] + 

- 15926803200 [Ax, [A 1; B 3 ]] - 3113510400^, [A l5 B 4 ]] - 359251200[Ai, [Ai,B 5 ]] + 

+ 821966745600A! - 92926310400B! - 503909683200B 2 - 366915225600B 3 - 142742476800B 4 + 

- 30656102400B 5 - 2874009600B 6 + 22800476160001 

(4.19) 

A 15 = 681080400H 15 + 20432412000H 13 + 286053768000Hn + 2479132656000H 9 + 

+ 14874795936000H 7 + 654491021 18400H 5 + 218163673728000H 3 - 6006[A n , [Ai, B x ]]+ 

- 2702700[A 9 , [Ax, BJ] - 540540[A 9 , [A 1; B 2 ]] - 665945280[A 7 , [A b BJ]+ 

- 242161920[A 7 , [Ai, B 2 ]] - 60540480[A 7 , [Ai, B 3 ]] - 84453969600 [A 5 , [A u Bi]]+ 

- 39956716800[A 5 , [Ai,B 2 ]] - 15437822400[A 5 , [Ai, B 3 ]] - 1816214400[A 5 , [Ai,B 4 ]]+ 

- 4544168428800 [A 3 , [Ai,Bi]] - 2332019289600[A 3 , [Ai,B 2 ]] - 1100625926400[A 3 , [Ai,B 3 ]]+ 

- 250637587200[A 3 , [Ai,B 4 ]] - 32691859200[A 3 , [Ai,B 5 ]] - 103981906828800[A 1; [Ai,Bi]]+ 

- 49669831411200^!, [Ai, B 2 ]] - 17980522560000[Ai, [A 1} B 3 ]] - 4729422297600[A X , [Ai, B 4 ]]+ 

- 893577484800[Ai, [Ai, B 5 ]] - 65383718400[Ai, [Ai, B 6 ]] + 560992303872000Ai+ 

- 5317875763200Bi - 354379753728000B 2 - 306780406732800B 3 - 148115916748800B 4 + 

- 44373750220800B 5 - 7758867916800B 6 - 610248038400B 7 + 17346736382976001 

(4.20) 

where we have used the following notation: 

2L 

H 3 = >/2^[e i ,[e i+1 ,e i+2 ]] (4.21) 
fc=i 

2L-2 

H 5 = V2 ^2 i e u N+i, [e< + 2, [ei+3, ei+4]]]] (4.22) 

k=l 
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2L-4 

H 7 = V2 ^ t e i> ( e i+i, [ei+2, [e i+ 3, [ei+4, [e i+5 , e i+6 ]]]]]] (4.23) 
fc=i 

B x = V2(e 2 + e 2L+2 ) (4.24) 

B 2 = V2(e 3 + e 2L+1 ) (4.25) 

B 3 = + e 2L ) (4.26) 

We call the A 2n _i the Lattice Local Integrals of Motion (IOM), since they play a completely 
analogous role as the I 2n _i of jl], as we shall realize even better in the next section. 
The A 2n -i form a maximal Abelian subalgebra (of hermitean operators) of the enveloping 
algebra of the TL algebra, in the same way as the CFT Local IOM I 2n _i span the maximal 
abelian subalgebra (of hermitean operators) of the enveloping algebra of the Virasoro alge- 
bra. As remarked in [9] the 2 enveloping algebras should become isomorphic in a suitably 
defined L — > 00 limit. 



5 Lattice Local Integrals of Motion 

In this section we want to obtain the eigenvalues of the Lattice Integrals of Motion which 
we have characterized in the previous section in terms of the Temperley-Lieb algebra. 
Starting from the following expansion: 

00 . 

logT(x) = + V , " (5.1) 
6 v ; ^ n\(2x) m v ; 

m=l v ' 

we readily obtain the following formulas for the eigenvalues of the A 2n _x: 

L+1 

A 2n = 2 2n (2n-l)\J2^ 2n (t k ) (5.2) 
fc=i 

L+1 

A 2n -i = 2 2n -\2n - 2)!(-2^sin 2 "" 1 (t fc ) + ^ sin 2 "- 1 ^)) (5.3) 

fceP k=i 

one then observes that the sum over sines in the previous equation can be evaluated analyt- 
ically thus yielding: 

X>-m =^B-irf; ') esc ( ( 2<„-*) - i )3I y. 

k=i k=o v 7 (5.4) 

/ / 4(n - k)L + 4(n - fc) + L 
• 1 — sin 7T 



4L + 6 

One then wants to expand the full eigenvalues A 2n _\ in powers of 4L + 6, for this reason one 
introduces the following coefficients: 

^ = ^j£ ( 2n 7 ') (-ir fe - 1 (2(n -k) - l) 2m - 1 (5.5) 

fc=0 ^ 7 



in terms of which one has the following expansions 

o n i, s o n ,(n-l)! 2 4L + 6 1 ^ 9m ,„ (-l) m / vr /l 

> sin 2 ™" 1 (t k ) = 2 2n ' i j J — + > 2 2m - 1 C n m - / — Vr( I B 2m { - 

^ y ' (2n-l) vr 2 ^ ' 2m !V4L + 6/ V2 

fc=l v y m=n v y 

(5.6) 

- 2 sin— (H) = £ £,„„ ^ (jJL_) " 4m(2* - 1)—' (5.7) 

m=n 1 

where B n (z) are the Bernoulli polynomials. Now, using this results one can obtain that 



kep k=i ^ >' 

m=n v y fceT 7 

in which we recognize that the following piece 

[AmJ2(k - -J 2 " 1 " 1 + = a m I 2m ^(V) (5.9) 

is proportional to the Local Integrals of Motion of the underlying CFT [B], where 

^n(2n - l)3 n r(4n - 1) , 

^2. ( — ]) m / 7T \2m-l 

A 2n _ x = + 2 2 "- 1 (2n - 2)! £ ^ m ^_2-(__j o^I^T?) (5.11) 

m=n 

where the divergent part is defined as: 

Af:_ x = 2^-i) ^- 1 ) !24L + 6 _ 2 2 (™- 1 )(2n - 2)! (5.12) 
[2n — 1) 7r 

notice that the first relevant term in A 2n -\ aside from the ones appearing in the divergent 
part is I 2n -i, thus motivating that the A 2n -i being called Lattice Integrals of Motion. 
We now remark that since C n ^ m = for m < n the m summation in the formula for A 2n -i 
can also be taken to start from m — 1. 

In order to obtain the explicit expression for the A 2n we now consider that: 

L+l 



(5.1 



k=l 

so that we have 



£^"(^-i(^ + i>-WV)) <" 3 > 



A 2n = -(2n - 1)! ( { ^(L + 1) - 2— + V ) ) (5.14) 



we now decompose the A 2n -i as: 

A 2n _ 1 =At v _ 1 + A 2n _ 1 (5.15) 

and resum the contributions of the divergent part of A2 n -i and of A 2n so that 

I TV \ _ ( \^ A 2n (L) \ - Afn-ljL) \ \ ^ ^2n-l fulfil 

° g W ~ (2x) 2 "(2ri)! ^ (2x) 2 ^ 1 (2n- 1)!/ ^ (2x) 2n ~ 1 (2n - 1)\ { ' 

n=l n=l v / \ ' n= l v / \ / 

where we have that 



£ (golfer ~ 1, !);(§,§ ); 1) - i log (i±f ) (5,8, 
furthermore the sum over the A 2n -\ can be rearranged as: 

£ (2x)^-H2n - 1)! = ^i^ilETsY a ^-i(r)P2 m -i(-) (5.19) 

n=l v / \ ' m=l v ' 

where 

m 2n-l 

-P2m-l(^) = -C nim (5.20) 



2n- 1 

n=l 



so that the expansion of the logarithm of T takes the form (analogous to the result of [7J for 
Polymers on the strip ): 



~ (-l) m / 7T \ 2m-l 1 

logT(x) = -2(2L + 3)f bulk (x) - f hou (x) + Y,7^S\{^T+3) a rnhm-i(V)P 2m -i(- 

m.=l ^ ' 

(5.21 

where we have introduced the bulk and boundary free energies: 



and 



= -j log (i(i + - i)) - ± 3*i((i i, i); (|, |); (5.22) 

= log (l - Jj) + i log (i±f ) (5.23) 

a; 

and by taking into account the relation: 

1 Is, 1 „ , / 1 .,3 3, 1 . 7T . / 1 . . 

log (^w + ;>* = x * FM r l ' *>< <? 2> ; ?) + 2 ^ + V 1 - ?> (5 ' 24) 

W) = ~ ( f log(-^y + §)* + log V2) (5.25) 

which is the well known formula for the bulk free energy of the Ising model. 

Notice also that the same kind of expansion has already been obtained for Critical Dense 

Polymers on the strip in [TJ. 



one gets 
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6 Characters and Operator Content 

Let p be the length of the partition V, one then defines the finitized characters to be 

x+(q)= Q^ er ^ (6.1) 



V, odd p 



X L M)= E (6-2) 
which, by defining the q-Binomial 



V, even p 



m—1 



m) J-l (i _ q i+i) 
7 1 i=0 v ^ y 

can also be expressed in the following Fermionic form: 



m— n N 



m=0 



L + l 

2m 



(6.4) 



X £( 8) = 5: " 2 "^ 2 "' + H 2 t- 1 i) (6 ' 5 » 

m=l ^ 7 9 

alternatively we also have the following Bosonic forms for the finitized characters: 

L+l L+l 



■A 



(9) = ^(11(1 + 9*-*) + 11(1 -A* - *)) (6.6) 
fe=i fc=i 

-, L+i l+i 
^-(9) = 2 ( + <? M ) " - <? M )) (6-7) 

k=l k=l 

We remark that these finitized characters have been already object of extensive study, and 
are found, in a slightly less symmetric form for example in pp. 

The characters with the + subscript enumerate the states corresponding to the b = +1 
boudary condition whereas the subscript — corresponds to b = —1. By use the above 
formulas one can see that the partition functions behave as: 

Z = TrfV 7 ) ~ Z div (L,M)q~^xi(q) (6.8) 

where 

_ M 

q = e 77 (2L+3)* (6.9) 

and: 

\ogZ dm = -2M(2L + 3)f bulk (x) - Mf bou (x) (6.10) 
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We now take a look at the expansions of the L = oo characters: 
x ~ ( q ) = 1 + q 2 + q 3 + 2g 4 + 2g 5 + 3g 6 + 3g 7 + 5g 8 + 5q 9 + 7q 10 + 8g n + 1 lq 12 + 12q 13 + ... (6.11) 

nI3579JAK|1517H)2123 

X- (^J = 9 2 +Q 2 +Q 2 + <7 2 +2g 2 +2g 2 +3g 2 +4g 2 +5g 2 +6g 2 +8q 2 +9g 2 + . . . (6.12) 
we notice that: 

x+(g) = xo(g) (6.13) 

X-(q) = Xi(q) (6-14) 

whre xo an d Xi are the respectively the characters of the irreducible h = 0, h = 1/2 
representations of the Virasoro algebra. 

The states belonging to the (0)+ and (~)_ representations, and labelled by a partition V, 
will all be of the form: 

\V) = ip gfep-i . . . ip 2fc 1 -i |o) (6.15) 
2 2 

with p respectively even or odd. 
The fermion field is expanded in this 

where the fermi modes satisfy in all cases the following anticommutation relations: 

{ipn^m} = Sn+mfi (6-17) 

All the above conformal characters and fermionic decompositions are well known and appear 
in many works, we point out here for example [5] pH] as particularly understandable. 
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